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Abstract 

We present numerical evidence for the spontaneous breaking of the Zf^ symmetry of four- 
dimensional twisted Eguchi-Kawai models with SU(A^) gauge group and symmetric twist, for 
sufficiently large N. We find that for N > 100 this occurs for a wide range of bare couplings. 
Moreover for < 144, where we have been able to perform detailed calculations, there is no 
window of couplings where the physically interesting confined and deconfined phases appear 
in the reduced model. We provide a possible interpretation for this in terms of generalised 
'fiuxon' configurations. We discuss the implications of our findings for the validity and utility 
of space-time reduced models as iV ^ oo. 



1 Introduction 



The conventional way of calculating the properties of the large N limit of SU(A^) gauge theories 
using lattice methods follows three main steps [IJ: the calculation of some dimensionless ratios 
of physical quantities at fixed lattice spacing a and fixed N, for several values of a and 
A^(= 2,3,4,5,...); the continuum extrapolation (a — > 0) of the resulting lattice values at 
fixed N] the large N extrapolation (A^ —>■ oo) of the resulting continuum values, assuming 
a leading correction of order (for pure gauge theories). This way of studying the physics 
of the SU(oo) theory has been the subject of intense research in the past decade P, and 
has allowed the accurate calculation of a number of nonperturbative properties of the pure 
gauge theory, such as the string tension pj, some glueball masses [3], and the deconfining 
temperature [3]. However, the numerical simulations involved become very costly when we 
consider gauge groups of size > 10 on reasonably large volumes. So, if one wishes to check 
by explicit calculation, at very large A^, that these extrapolations are indeed under control, 
one needs another approach. Fortunately such an alternative approach exists: it uses an old 
idea introduced by Eguchi and Kawai that goes by the name of reduced models O [6] . 

In their original proposal [5j, Eguchi and Kawai showed that an exact correspondence can 
be formally constructed between the N = oo limit of SU(A^) gauge theories living on an infinite 
(i-dimensional lattice (the original theory) and the same theory living on a 1*^ periodic lattice 
(its reduced model). The correspondence states that Wilson loops in the original theory and 
its reduced model obey the same Dyson- Schwinger equations, and therefore it is reasonable to 
deduce that their numerical values should coincide (ignoring the possibilities raised by multiple 
solutions). However, these loop equations are only identical if we make two assumptions: 1) 
the large N factorisation of the physical observables, and 2) an unbroken symmetry of the 
reduced model. While it is reasonable to believe that factorisation is a universal property of a 
large class of physical operators in the large limit [7], nothing prevents the symmetry 
of the reduced model from being spontaneously broken. Indeed, the symmetry was found to 
be spontaneously broken in the continuum limit of the original Eguchi-Kawai model [8J, so 
the Eguchi-Kawai correspondence is spoiled in that case. But there are alternative ways of 
constructing reduced models of SU(A^) gauge theories that keep the Zf^ symmetry unbroken 
for arbitrarily weak couplings: 1) by quenching the eigenvalues of the reduced link matrices 
(quenched Eguchi-Kawai model, or QEK model) [$], or 2) by imposing twisted boundary 
conditions on the reduced lattice (twisted Eguchi-Kawai model, or TEK model) [9]. 

A great amount of work was done in the first half of the 1980s in trying to extract phys- 
ical properties of large A^ theories from numerical simulations of (mainly) TEK models [B]. 
A drawback of these models, however, is that the volume of the lattice on which the original 
theory is defined is directly related to the number of internal degrees of freedom of the reduced 
model, so very large gauge groups must be used in order to obtain information that is suffi- 
ciently near the continuum and large volume limits to be physically relevant. Consequently, 
due to obvious computational constraints, the results obtained in the early days lacked ac- 
curacy and, in general, did not allow one to reach reliable conclusions about the physics of 
the N = oo theory (for example, we are not aware of any finite-temperature study using this 
framework that demonstrated, in a single calculation, a clear separation between the bulk 
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and deconfining phase transitions). With present day computational resources, one can easily 
perform more precise calculations and for larger gauge groups than before, and this is one 
motivation for returning to a careful study of large N physics using reduced models. There 
are other motivations for such a study: first, it allows us to do simulations with much larger N 
than with the traditional method, since the number of external degrees of freedom (volume) 
is reduced to its minimum; and second, it provides an alternative method to approach large 
physics, with which we can compare (and hopefully confirm) the numerous results already 
obtained via the conventional method. Finally we point to the importance of reduced models, 
in particular the twisted version, in matrix models of noncommutative field theory [10] and 
M-theory [11]. 

In this paper, we will summarise the results of a detailed numerical study of the properties 
of symmetric-twist SU(A^) TEK models in four dimensions. In the next section, we provide a 
short review of TEK models. Next, we present the details of our numerical simulations and 
discuss the results obtained from them. The surprise is that we find that the Zf^ symme- 
try of the TEK model, usually assumed to be unbroken for all couplings, is in fact broken 
spontaneously at intermediate couplings when is sufficiently large. At the same time there 
appears to be a persistent metastability at weak coupling that prevents the appearance of a 
physically interesting phase in the reduced model. We conclude with a discussion about the 
possible origin and the implications of this symmetry breaking. Some of these results have 
been presented at SMFT 2006 [12]. 



2 TEK model 

The SU(A^) TEK model in four dimensions is defined by the partition function on a periodic 
lattice with a single site, 

Ztek = [ [dU] e-^^^™- (1) 



where [dU] = Y[^=i ^^f^ a product of SU(A^) Haar measures, t/^ G SU(A^) are the link 
variables on this reduced lattice, b = l/g'^N is the inverse bare 't Hooft coupling and S'tek is 
the action of the TEK model, 

4 

Stek = STEKin; [U]) =J2 (/ - z^,{n)U^UMlUl) > 0. (2) 

S'tek depends on the reduced link variables and on the twist-tensor, n^,^ = ~n^^ ^ '^n, 
that labels the twist factor z^i,{n) = e'^^^'^" G Z^. The twist-tensor defines uniquely (modulo 
SL{4:, Z) transformations) the geometry of the effective lattice on which the original theory is 
defined. We choose to use the standard symmetric twist, n^i. = V^, V /x > i/, of Gonzalez- 
Arroyo and Okawa [9] . With this twist the TEK model corresponds to an effective symmetric 
L"^ lattice with A^ = L^. Under certain conditions, as discussed below, the TEK model has the 
same planar limit as a conventional SU(A^) lattice theory with a plaquette action, at the same 
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value of the 't Hooft bare coupling. This restricts the physical value of to be the square of 
an integer. 

The lattice operators on the reduced 1^ lattice are obtained by replacing the link variables 
in their original definition with the corresponding reduced link variables, and then introducing 
the appropriate twists. The I x J Wilson loop in the {/i, u} plane in the TEK model is then 
given by: 

W,M, J) = ^.uinY'UlU^Uj/uy (3) 
and the Polyakov loop in the /i direction (on the equivalent L = \/N lattice) is given by: 

P, = Ui (4) 

The action of the TEK model in ([2]) has two distinct symmetries acting on the link vari- 
ables: the usual gauge symmetry, 

^ nUf,n\ n e sv{n) (5) 

and an additional Zjj symmetry, 

Uf, ^ z^Uf,, z^ e Zn (6) 

In order for the TEK model to be able to reproduce the same large physics as the original 
theory, a necessary condition is to have vanishing expectation values for the traces of reduced 
open lines, such as f/^ with a not an integer multiple of L = In the original theory this 
is guaranteed by the fact that such an open line is not gauge invariant, but in the single-site 
reduced model any such operator is trivially gauge-invariant and could be non-zero. However, 
unlike Wilson loops ^ and Polyakov loops (jl]), we note that reduced open lines are not 
invariant under the Z^ symmetry (jS]). So, if the Zj^ symmetry is not spontaneously broken, 
the expectation value of open lines will indeed be zero, and the loop equations of the TEK 
and conventional theories will be identical in the planar limit. 

In the strong coupling limit, 6 — >■ 0, the partition function in eqn([T]) is dominated by the 
Haar measure [dU]\ consequently, the expectation value of the trace of a link variable (and in 
general of any open line) is automatically zero, so the model has an unbroken symmetry in 
that limit. We call this phase the random phase and denote it by Z^j^'^\ In the weak coupling 
limit, h —>■ oo, the partition function in eqn([T]) is dominated by the classical vacuum of the 
TEK model, i.e. by the field configuration f/^ = that solves the equation: 

z,,in)T,T,TlTl = I. (7) 

The matrices defined above, also known as twist- eaters, saturate the lower bound of the 
TEK action in eqn([2]). They have the property of being traceless, as do most of its powers: 
Tr(r") = 5o,(a mod L)- In other words, the classical vacuum of the TEK model has a Zfj- sym- 
metry that ensures that all lines (except Polyakov loops) vanish. Small quantum fluctuations 
around the twist-eating configuration also keep the Zff symmetry unbroken. We call this phase 
the twist-eater phase and denote it by Zfj-^\ Extrapolating from the fact that the symmetry 
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is unbroken in the 6 — *• and b —>■ oo limits, it is generally claimed that the Zf^ symmetry 
of the TEK model should be unbroken for all couplings. There is, however, no compelling 
theoretical argument justifying such a claim for the intermediate-coupling regime, although it 
receives some support from the fact that previous numerical simulations of the model [6] have 
not shown any sign of symmetry breaking. 

3 Results 

In this section we provide numerical evidence that for sufficiently large A^, namely > 100, 
the Zff symmetry of the TEK model is spontaneously broken. For iV < 81 we see no sign of 
the symmetry breaking, which is consistent with the results from numerical simulations done 
in the past [6]. 

3.1 Phase structure 

In the large limit, SU(A^) lattice gauge theories possess a set of phase transitions that we 
would expect to observe in our simulations of the TEK model, if the Eguchi-Kawai correspon- 
dence between the planar limits of the reduced and parent lattice gauge theories is in fact 
valid. First, as we increase b we expect to encounter a first-order phase transition between the 
strong and weak coupling phases of the theory. If the lattice is large enough the transition is 
between two confining phases. As we increase b further, the lattice will encounter a deconfining 
phase transition. This will occur when the smallest lattice length satisfies a{b)L^ = 1/Tc, 
and it is characterised by the Polyakov loop in that direction acquiring a non-zero vacuum 
expectation value As we increase b further there will be a sequence of three further transi- 
tions corresponding to each of the other three Polyakov loops acquiring a non-zero expectation 
value. Finally we are in a phase where the volume is small and the physics on all length scales 
is perturbative and there are no further transitions as a — > 0. Although it is not obvious that 
it should be so, it has been observed [H] that exactly the same phase structure persists on a 
symmetric lattice, and hence should be observed in the TEK model. We now discuss these 
transitions in more detail. 

The bulk transition is a lattice artifact that occurs, for the standard plaquette action, 
a.t b = bs ^ 0.36 in the large A^ limit. It separates the unphysical phase at strong lattice 
coupling (known as the bulk phase) from the physical phase at weak lattice coupling, through 
which one approaches the continuum limit of the theory. The bulk transition also occurs in 
the TEK model, in this case separating the random phase, Z^^\ from the physically relevant 
twist-eater phase, Z^^^; it was accurately determined in [T5| to have the value bs = 0.3596(2) 
in the N = oo limit. 

For large enough lattices, the deconfining transition is located on the physical side of the 
bulk transition, separating the low T confined phase (with non-zero string tension, a ^ 0) 
from the high T deconfined phase (with a = 0). It is signaled by the spontaneous breaking 
of the Zn symmetry of the Polyakov loop along the direction of the 4-torus with the 
shortest physical length, aL = T~^, when the temperature T rises above a critical value 
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Tc] it is also accompanied by a small discontinuity in the average value of the plaquette 
(reflecting the latent heat of this first-order transition [^). The critical coupling he at which 
the deconfining transition occurs scales with the lattice size L, so it is natural to expect that 
there is a critical size L = Lb below which the deconfining transition is overtaken by the bulk 
transition, and only for L > Lb are the two transitions decoupled. From the calculations of 
Tc versus N in [4] we can infer that ~ 9. In the TEK model, we can therefore expect 
to observe a corresponding deconfining transition, signaled by the spontaneous breaking of 
the Zb symmetry of the reduced Polyakov loops [17], once N = L^ > ~ 81. There 
is some uncertainty in this value, that will be enhanced in practice by any metastabilities 
associated with these first-order transitions, so it is only at very large that we can expect 
the symmetric-twist TEK model to exhibit a genuine weak coupling confined phase. 

The TEK model is equivalent to a symmetric lattice, with L = \/N, and so there is no 
direction that naturally plays the role of an inverse temperature. This system has been studied 
in [131 . One finds that as one increases b from the confining phase one encounters a sequence 
of four transitions at each of which Polyakov loops along one more direction acquire non-zero 
expectation values. The first of these transitions occurs when a{b)L = and is indeed just 
the deconfining transition (as follows from the observed vanishing of finite volume corrections 
when N ^ oo The phases are labelled in [13] by Xc {X = 0, . . . , 4), with Oc the confining 
phase and 4c the small volume phase. Each of of these transitions has a continuum limit 
and hence a specific physical interpretation. Therefore, if the Eguchi-Kawai correspondence 
is correct, we should also be able to observe these Xc phases in the symmetric-twist TEK 
model, on the physical weak coupling side of the bulk transition. 

3.2 Simulation details 

We have performed Monte Carlo simulations of the four- dimensional SU(A^) TEK model 
with the conventional symmetric twist tensor (n^^, = L = a/ZV, for fx > u). The thermalised 
configurations at fixed b were generated using the heatbath algorithm constructed by Fabricius 
and Haan [18]. As a check, a subset of calculations were performed with a straightforward 
(albeit less efficient) Metropolis algorithm, specifically a modified version of the algorithm 
originally constructed by Okawa for the untwisted model |19]. The simulations started either 
from a classical twist-eater configuration, referred to as a cold start, or from a completely 
randomised configuration (in practice one that was thermalised at 6 = 10~^) referred to as a 
hot start. 

We considered gauge groups = = 25, 36, ... , 144 corresponding to effective lattice 
sizes between L = 5 and L = 12. For each case, we performed two different simulations differing 
in the range of values of b spanned: in the first run (run A), the interval b G [0.10,2.00] was 
spanned with large steps (varying between |A6| = 0.01 and |A6| = 0.02); in the second run 
(run B) we wanted to analyse in more detail the region where the phase transitions occur, 
so we spanned the interval b G [0.200,0.500] with smaller steps, |A6| = 0.005. The other 
parameters of the simulations, like the number of sweeps or the number of measurements, 
varied with each case and are summarised in Table [1] for the larger values of A^. 

We calculated the expectation values of several reduced observables, namely the traces and 
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eigenvalue densities of Wilson loops of several sizes, of Polyakov loops, of link variables (f/^) 
and of some other open lines {U^, for a = 2, . . . , 5). 

3.3 Results for iV < 81 

The phase structure of the TEK model for < 81 is very simple. Essentially, there is only 
one phase transition, at 6 ~ 0.36 [15\. The transition is strongly first-order, as demonstrated 
by its relatively large hysteresis cycle, and in each one of the two phases the average plaquette 
fits very well the strong- and weak-coupling expansions of both the TEK model and Wilson's 
lattice gauge theory, as shown in Fig. [H Moreover the location of the transition is consistent 
with that of the bulk transition in usual lattice calculations at large N. So we can confidently 
identify it as being the bulk transition. 

We see in Fig. |2] that the average traces of the link variables vanish in the whole range 
of lattice couplings (although the fiuctuations at intermediate couplings are relatively large). 
Therefore, we can conclude that the Eguchi-Kawai correspondence between the planar limits 
of the reduced model and the parent (Wilson) lattice gauge theory is valid at least up to 
= 81. As an example of this, we show in Fig. [T]not only the 3-loop perturbative result for 
the plaquette in Wilson lattice gauge theory at A^ = oo, given by [T6] : 

, , 1 0.653687 0.4066406 

but also some values that we obtained from a direct simulation of Wilson lattice gauge theory 
for 4 < A^ < 16, extrapolated to A^ = oo. We see from Fig. [1] that these extrapolated 
values coincide with the TEK values, within small errors that are consistent with the various 
corrections and systematic errors. 

The phases separated by the bulk transition can be identified with the random phase, Z'l^'^\ 
and the twist-eater phase, Zf^^\ The justification comes from an analysis of the eigenvalue 
densities of the link variables: in the random phase, the eigenvalues of the link variables are 
uniformly distributed over the unit circ 1^ , while in the twist-eater phase the eigenvalue density 
of link variables is composed of L = separated lumps spread around the elements of Zl 
on the unit circle (modulo a global Zn rotation of the whole spectrum). The latter is due to 
fiuctuations around a twist-eating configuration (modulo a Zn rotation), whose eigenvalues 
are the elements of Zl [20\ . 

However, for A^ < 81, it is not possible to pinpoint a deconfining transition in the twist- 
eater phase, since the Polyakov loop and the plaquette suffer a discontinuity only at the bulk 
transition. This result should not be too surprising, since, as remarked earlier, the sizes of the 
effective lattices involved, L < 9, are probably smaller than the critical size Lb needed to have a 
clear separation between the deconfining and bulk transitions. Since Polyakov loop expectation 
values are all non-zero in the weak coupling twist-eater phase, while the strong coupling phase 
is confining, the bulk transition can easily be misinterpreted as being a deconfining transition, 
particularly in calculations with low statistical accuracy. The only way to be sure that one 



^Due to unitarity, the eigenvalues of U(iV) or SU(-/V) matrices are constrained to be elements of U(l). 
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has not confused the two transitions is to vary b so as to identify both transitions exphcitly 
in a single calculation. As we have remarked earlier, none of the old TEK calculations that 
claimed to find a deconfining transition, displayed such a result. We are therefore forced to 
discount those early claims. 

This negative result motivated us to extend our investigation to L > 10, whose results are 
summarised next. 

3.4 Results for N > 100 

The phase structure of the TEK model for > 100 turns out to be much richer than for 

< 81, involving several transitions to new phases. The main characteristic common to 
all those new phases is that the Zfj symmetry is spontaneously broken, thus violating the 
Eguchi-Kawai correspondence. 

In the A^ = 100 case, the average value of the real part of the plaquette undergoes two 
distinct transitions as we decrease b from a cold-start, but only one transition occurs when 
we increase b from a hot-start, as shown in Fig. [3l In the case of the cold-start scan, the 
phase in between the two transitions is characterised by a non-zero expectation value for the 
trace of one of the link variables as shown in Fig. H] (and also by an imaginary part of the 
plaquettes containing that link variable). This implies that the Zff symmetry of the TEK 
model is spontaneously broken to a Zf^ symmetry in the intermediate-coupling phase (the 
traces of the link variables still have zero expectation value in the other three directions). As 
discussed at the end of section [21 this effect should not come as a complete surprise, since there 
are no compelling theoretical arguments that prohibit the breaking of the Zj^ symmetry at 
intermediate coupling^. The transitions to the regime with Zff symmetry have a discontinuous 
character, indicating that these are first-order phase transitions (or would be if the number of 
degrees of freedom were infinite). Increasing b from a hot-start only one transition is observed, 
and that is the one from the random phase to the Z^-symmetric phase. As we see in Fig. |5l 
there is no transition from the Zf^ phase to the expected weak coupling twist-eater phase, even 
at very large values of the inverse lattice coupling. This implies the existence of stable extrema 
of the TEK action that are Zff non-preserving and from which we are unable to tunnel because 
of large barriers. This immediately raises the question whether this effect persists for larger 
A^ and what consequences it has for the Eguchi-Kawai correspondence. 

As at lower A^, there is no sign of a deconfining transition in the physically relevant phase 
of the SU(IOO) TEK model. By 'physically relevant' we mean the twist-eater phase, since it is 
the only phase that we observe having the necessary Zf^ symmetry while being continuously 
connected to the continuum limit. 

For A^ = 121 and A^ = 144 we find, at intermediate lattice couplings, a sequence of transi- 
tions, as illustrated in Fig. [6] and in more detail in Fig. [71 These transitions are similar to the 
one we saw for A^ = 100, in that each is associated with the breaking or restoration of one or 

^That the imaginary part of the plaquettes can be non-zero follows from the fact that the action of the 
TEK model in ([2]) is not CP-invariant. At small b it can be calculated in a strong coupling expansion and one 
finds that the contribution decreases with N, so that it is real in the planar limit and is consistent with the 
Eguchi-Kawai correspondence. 
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more of the symmetries. Thus they lead to non-zero expectation values for the correspond- 
ing link variables, as shown in FiglHl We can label the phases by the number k of preserved 
Zjv symmetries, as (/c = 0, . . . , 3). In the = 144 case, all transitions between different 

phases are clearly separated, while for = 121 the simultaneous breaking/restoration of 
several Zn symmetries can occur. Tables [2] and [3] summarise the critical couplings at which 
the observed transitions occur. 

Just as for N = 100, when we increase b from strong coupling for the SU(121) and SU(144) 
TEK models, we observe no tunnelling from the final Z^ phase to a twist-eater phase. This 
reinforces the idea that some stable extrema of the TEK action are contributing to the breaking 
of the Z^ symmetry. The nature of these extrema can be inferred from the eigenvalue spectrum 
of the link variables. In a phase, the k link variables of the Z^r-preserved directions have 
a uniform eigenvalue density, while in the other N — k broken directions their eigenvalue 
density consists of a bell-shaped lump with finite support, spread around an element of Z^, 
as illustrated in Fig. O At larger b the lump becomes narrower, as shown in Fig. [TDl indicating 
that in the limit 6 ^ oo all eigenvalues collapse to one value which is an element of Z^y. In 
sum, the numerical simulations suggest that the breaking of the Z^ symmetry in the TEK 
model for A^ > 100 might be caused by fluctuations around (stable) configurations of the form 



Our numerical calculations suggest that the gauge fields responsible for the spontaneous break- 
ing of the Zfj- symmetry are of the form of fluctuations around G Z^r. These field config- 
urations appear to be able to survive reasonably large fluctuations at intermediate couplings, 
and to become increasingly stable as we increase b to arbitrarily large values. This requires 
the existence of stable extrema of the TEK action that break the Z^ symmetry. At the same 
time we see explicitly from Figs. [3] and E] that their action is substantially larger than that of 
the 'twist-eating' fields obtained at the same b from a cold start. Thus the (infinite?) metasta- 
bility requires large barriers, which implies the existence of corresponding unstable extrema 
of Stek- We will now argue that the elements of Zn may provide such configurations. Our 
arguments below will follow very closely the calculations of van Baal in his study of surviving 
extrema of the TEK model [21 j. 

Consider a small fluctuation around a general configuration Q^, parameterised by an (anti- 
hermitian) element of the algebra of SU(A^): 



The expansion of the TEK action ([2]) around fl^ up to second-order in is given by: 



e Zn. 



4 Interpretation 



(9) 
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where P^^, = z^^{n)Qjj^i}lflfj^fl^ is the reduced plaquette, P{ij,u} = | {Pfiu + Pu^i) is its symmetri- 
sation, P[^,^] = ^ (P^j^ — Pj,^) is its antisymmetrisation, D^X^ = Q'^^X^Q^ — Xy is a discretised 
covariant derivative, and P^jy = D^Xy — D^X^j, is a discretised field strength tensor. The sta- 
tionarity condition for the configuration Vt^ to be an extremum, ^S'tek = 0, then imphes the 
equation: 

4 4 

<55tek = - I] Tr (P[H^M-) =0^5^ (^^.P[H^l " ^[h) =0 (H) 

An obvious class of solutions to the stationarity condition above are the configurations sat- 
isfying Pfj^y G Z]^. In fact, van Baal gave strong arguments [21] for the conjecture that all 
stable extrema of the TEK action have this form. So let P^v = z^y{n — m), where n = n^j, is 
the twist-tensor and m = rrif^u is another integer-valued matrix. Expanding the TEK action 
around a configuration from this class reduces (|T0|) to: 

Tr(-P^2j + 0(X3) (12) 

The stability condition, 6'^Stek > 0, is equivalent to: 

> 0, V/i,z/ (13) 

since Tr(— P^^) = O(X^) is always non-negative. 

Let us consider particular solutions to the stationarity condition P^^ = z^^^n — m), i.e. 
Zfj,,^{Tn)fljj^Qlfl^Qi, = I. For the case m^j, = n^^, the unique solution to P^^, = I is obviously 
the twist-eating configuration, the classical vacuum of the TEK model. Now consider the case 
when m^i, ~ 0{^/N) and (fT3|l is satisfied. In this case the solutions to P^^, = Zf^i^{n — m) 
(if they exist) are called fluxons [2T]. Fluxons survive the large limit, being stable minima 
of Stek, and share some properties with twist-eaters, like tracelessness; therefore, fiuxons 
are not expected to contribute to the breaking of the Zf^ symmetry. Next consider the case 
rrifj^iy = 0. The solution to P^u = z^^^n) is the set of all SU(iV) diagonal matrices. In general, 
these matrices have non-zero trace, and so they are good candidates for being the stable 
extrema that break the Z'^ symmetry. As a consequence of the first term on the l.h.s. of 
( IT2|) . the action of the diagonal matrices grows proportionally with A^, i.e. in the partition 
function such diagonal matrices appear to be suppressed relative to twist eaters by a factor of 
Q-bN ^ While naively this appears to be a large suppression, in a system with 0{N'^) degrees of 
freedom this is not so. Quantum fiuctuations will naturally contribute terms of 0{N'^) to the 
effective action which can easily overwhelm the difference in the classical action. At fixed A^, 
for large enough 6, the action difference eventually dominates; but if the barrier between these 
diagonal fields and the twist-eating fields has a natural 0{N^) value, then the tunnelling to 
the twist-eater phase would be suppressed by a factor like e~^^^ and would essentially never 
occur in a finite simulation - just as we appear to see in Figs. [3l [6] and [3 Such a barrier 
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will be related to a maximum or saddle-point of the TEK action whose action grows oc N'^. 
Possible candidates for such maxima would be configurations satisfying P^^ = z^^{n — m), 
with m^jy ~ 0{N) in such a way that the condition ( fT3l) is violated. 

At the classical level, the set of all diagonal SU(iV) matrices form a degenerate set of 
extrema of the TEK action. However, if we consider the effect of quantum fluctuations, we can 
easily see that the degeneracy is lifted and only the elements of Z^- become the true extrema. 
This is so because, for arbitrary X^, D^X^ = for elements of while D^X^ ^ for other 
diagonal matrices. Consequently the O(X^) term in (fT2l) . which is oc Tr(— F^^), vanishes in 
the former case and is always positive in the latter. Therefore, when quantum fluctuations 

are taken into account, the action of a general diagonal matrix becomes larger than the 
action of an element of Zjv, if h is large enough for 0{X^) -C O(X^), and the degeneracy 
is lifted. Therefore, the elements of provide stable vacua that break the symmetry. 
This mechanism is in fact the same as that which leads to the Zj^ symmetry breaking in the 
original untwisted Eguchi-Kawai model [22]; it is also similar to mechanisms that generate 
'order from disorder' in condensed matter physics. 

Let us therefore compare the contributions to the TEK action of the fluctuations X^ 
around twist-eaters (F^) and around elements of Z^ (z^), 



We might speculate that the different sensitivities with respect to the fluctuations might 
be at the origin of the spontaneous breaking of the Zff symmetry at intermediate couplings. 
In this regime the fluctuations are reasonably large, but the Haar measure does not dominate 
yet. The contribution to the action of the fluctuations around twist-eaters is larger (~ 0{X^)) 
than for elements of Z^^ (~ 0{X^)), so it might be possible that at a certain critical value of 
the coupling b the fluctuations around elements of Z^v will cost less action than the fluctuations 
around twist-eaters, and consequently the elements of Zj\f become preferred over the twist- 
eaters. This argument is, however, very speculative, since it uses perturbation theory to justify 
effects that are observed at intermediate-to-strong couplings, where phase transitions are also 
involved. 

Returning to the sequence of transitions observed in the large N TEK model, as shown in 
Fig. [7] for N = 144, we note their qualitative resemblance to the transitions described in [H] 
that occur on conventional symmetric lattices for large X as 6 — oo. This might suggest a nat- 
ural identification between the Xc phases and our Z^ phases, namely Xc = Z'^^ . However, 
there are at least two fundamental differences that undermine such an identification. First, 
the reduced model version of the so-called Oc phase (the confined phase on the physical, weak 
coupling side of the bulk transition, which is situated between the strong coupling random 
phase and the Ic phase where Polyakov loops in one direction have a non-zero expectation 
value) was never observed in any of the simulations. Second, and more importantly, our tran- 
sitions do not seem to scale with L = a/X in the same way as do the ones in lUj. That is 
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(14) 



(15) 
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to say, the critical couplings of the latter increase with the lattice size L in such a way that 
the Xc phases have a proper continuum limit. However, according to Tables [2] and [3l our 
transitions do not show any sign of the scaling with N that would be expected if we were 
seeing a physical transition on the equivalent L = \fN lattice that occurred at a fixed value 
of the physical length I = aL = a^/N. That is to say, they appear to be lattice transitions 
rather than physical transitions. Moreover all our phases are situated below the value of 
b ~ 0.36 expected for the bulk transition, and not above (as in the case of the Xc phases), and 
finally their breaking of the centre symmetry undermines any connection with a supposedly 
equivalent conventional lattice. On the other hand, if we forget the Eguchi-Kawai correspon- 
dence (which, after all, we have shown is invalid in this range of b) and consider the TEK 
model simply as a l'^ lattice, then the links are Polyakov loops, and the identification with 
the Xc phases of |13] is essentially trivial. This suggests looking at the untwisted 1^ lattice, 
where such phases are known to occur [H]. In Fig{TT]we show a plot of this phase structure 
that we have calculated for = 81. The phase structure we see is essentially identical to that 
obtained in the large TEK model with a hot start, as shown in Fig. [TJ and the interpreta- 
tion of the phases appears to be identical. It appears that for large enough A^, once quantum 
fluctuations are sufficiently important, the effect of the twist factor in the symmetric TEK 
model, z^jj{n) = e^'^*/^, vanishes for all practical purposes. 

5 Conclusions 

We have analysed numerically the properties of four- dimensional SU(A^) TEK models with 
a conventional symmetric twist [9]. The gauge groups used in our simulations ranged from 
= 25 to = 144. 

For A^ < 81, the properties of the TEK model are essentially the same as those obtained 
in earlier simulations. As we vary the inverse bare 't Hooft coupling, b = 1/g'^N, we encounter 
a single first-order phase transition separating the strong and weak coupling regions and, 
crucially, the centre symmetry which is needed for the Eguchi-Kawai correspondence to be 
valid, is indeed maintained at all b. However our greater accuracy also makes it clear that there 
is no confining phase on the physical weak coupling side of the bulk transition, in contrast 
to some claims in the early literature. We pointed out that this should not occasion great 
surprise since we know, from conventional calculations of the deconfining temperature in 
SU(A^) gauge theories [4], that the length of the equivalent lattice, I = aL = a^/N, is probably 
less than 1/Tc for all b in the weak coupling region. 

For A^ > 100, on the other hand, we observe a clear breaking of the Zfj- symmetry over a 
wide range of intermediate couplings, challenging the conventional wisdom that this symmetry 
is unbroken for all couplings. We observe that the breaking of the symmetry occurs step 
by step in a sequence of transitions, with intermediate phases in which one or more of the Zn 
symmetries are unbroken. When, for our largest values of A^, we increase b from 6 ~ the final 
weak coupling phase appears to be one in which all four centre symmetries are spontaneously 
broken. At large enough b it becomes clear that the symmetry breaking fields are nothing but 
fiuctuations around centre elements. These are stable minima of the twisted lattice action, as 
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originally discussed by van Baal [21]. Since the action of these 'generalised fluxons' is only 
0{N) higher than that of the absolute minimum (the twist eater), it is not surprising that 
it does not tunnel to the latter since the generic barrier height will be O(iV^) in an SU(iV) 
gauge theory. 

The sequence of symmetry breaking transitions in the TEK model at very large is 
essentially identical to that which occurs in the original untwisted Eguchi-Kawai model (apart 
from a shift in the location of the transitions). It is as if the effect of the twist, a phase factor 
^2m/^/N associated with each plaquette, becomes irrelevant for large enough A^. Naively one 
might expect this to happen once the phase is small, i.e. for l-n j ^J~N <^ 1 or, equivalently, 

^ 477^. This provides a very large scale for what what we mean by being 'large', and 
might explain why new effects arise for A^ > 100 in the TEK model, in contrast to conventional 
calculations where finite A^ corrections are, typically, already very small for A^ = 6 or A^ = 8. 

One might still hope to encounter physically interesting phases by beginning at 6 = oo 
with a twist-eating field configuration and then decreasing h. If the first transition to strong 
coupling were to be at 6 ~ 0.360, as it is for A^ < 81, then the 12^ lattice that is equivalent to 
the SU(144) TEK model, should be large enough to be in the physical confining phase before 
reaching that value of h. Frustratingly, however, for SU(144) the critical value of h increases, 
as we see in Tables |2] and [3l so that there is once again no window for any weak coupling 
confining physics. Whether this will continue if A^ is increased further is obviously a very 
interesting question. 

On conventional symmetric lattices at large A^ there are additional physical phases 
for aL < 1/Tc, corresponding to the Polyakov loops acquiring non-zero expectation values 
one after the other ^14J. The first of these, as one increase b, is simply the usual deconfining 
transition and will occur at aL = 1/Tc. The subsequent transitions are expected to occur at 
higher b [H] although there are unfortunately no values in the literature for the corresponding 
critical values of b. It is not clear why we have not observed any of these transitions when we 
reduce b from the twist-eating start, particularly in the A^ < 81 case where the Eguchi-Kawai 
correspondence should hold for all b. It may be that there is a large metastability so that the 
twist-eating configuration tunnels to the bulk phase (for A^ < 81) or to the symmetry broken 
phase (for larger A^) before it can tunnel to one of these physical 'small volume' phases. Since 
the Eguchi-Kawai equivalence only holds for the planar limit, one cannot expect it to extend 
to the details of transitions and metastabilities which are typically driven by effects that are 
exponentially small in A^, and so there may be some surprises here. 

In this paper we have focused on two particularly interesting results from our study of 
the twisted Eguchi-Kawai model: the sjTiimetry breaking that invalidates the Eguchi-Kawai 
planar correspondence over a wide range of couplings, and the apparent inaccessibility of any 
physical phase even at A^ = 144. We have been able to make substantial analytic progress in 
understanding these phenomena, using configurations that are a generalisation of the 'fluxons' 
studied long ago by van Baal. Even if from a practical point of view our results undermine the 
utility of the most conventional TEK models, it would be theoretically interesting to pursue 
the numerical calculations to somewhat larger A^, so as to obtain some insight into what 
happens at asymptotic A^, and to develop further our analytic understanding of the model. 
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This we hope to do elsewhere where we will also present our results for TEK models with 
anisotropic lattice spacings and for the partially reduced models that one would need to use 
for calculations of the planar mass spectrum. 

Note added: During the course of this work we became aware in [23] of an unpublished 
talk in which observations were made about centre symmetry breaking at intermediate 
couplings in TEK models. 
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Table 1: Parameters of the numerical simulations performed at each value of the inverse 't 
Hooft coupling, 6; N is the size of the gauge group; nheats is the number of thermalising sweeps; 
Usweeps IS the number of sweeps used for measurements. 
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Table 2: Critical values of the inverse 't Hooft coupling, 6, associated with the breaking of one 
(or more) symmetries of the TEK model. This table refers to simulations that started from 
randomised configurations. refers to a phase of the TEK model with k unbroken directions, 
while the arrows mean transitions between those phases; a numerical value spanning multiple 
columns corresponds to the critical coupling associated with the simultaneous breaking of 
more than one Zjv symmetry. 
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Table 3: Critical values of the inverse 't Hooft coupling, 6, associated with the breaking (or 
restoration) of one (or more) Z^ symmetries of the TEK model. This table refers to simulations 
that started from twist-eating configurations. Z% refers to a phase of the TEK model with 
k unbroken directions, while the arrows mean transitions between those phases; a numerical 
value spanning multiple columns corresponds to the critical coupling associated with the 
simultaneous breaking (or restoration) of more than one Zj^ symmetry. 
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Figure 1: Average value of the real part of the plaquette, (Re Up), in the SU(64) TEK model 
versus the inverse bare 't Hooft couphng, b. The squares (□) represent the — > oo extrap- 
olation of the average plaquette in Wilson's lattice gauge theory. The solid line ( — ) and the 

dashed line ( ) represent the strong-coupling and the 3- loop weak-coupling expansions of 

Wilson's lattice gauge theory, respectively. 
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Figure 2: Average value of the real and imaginary parts of traced link variables, -^Tr f/^, in the 
SU(64) TEK model versus the inverse bare 't Hooft coupling, b (from a cold start simulation). 
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Figure 3: Average value of the real part of the plaquette, (Re -Up) in the SU(IOO) TEK model 
versus the inverse bare 't Hooft coupling, h. The squares (□) represent the N ^ oo extrapo- 
lation of the average plaquette in Wilson's lattice gauge theory. 
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Figure 4: Average value of the real and imaginary parts of traced link variables, j^Tr Un, 
in the SU(IOO) TEK model versus the inverse bare 't Hooft coupling, b (from a cold start 
simulation); the dot-dashed line (- • -) represents the average real plaquette. 
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Figure 5: Average value of the real and imaginary parts of traced link variables, ^^Tr U^, in the 
SU(IOO) TEK model versus the inverse bare 't Hooft coupling, h (from a hot start simulation); 
the dot-dashed line (- • -) represents the average real plaquette. 
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Figure 6: Average value of the real part of the plaquette, (Re -Up), in the SU(144) TEK 
model versus the inverse bare 't Hooft coupling, h. The squares (□) represent the N ^ oo 
extrapolation of the average plaquette in Wilson's lattice gauge theory. The solid line ( — ) and 

the dashed line ( ) represent the strong-coupling and the 3-loop weak-coupling expansions 

of Wilson's lattice gauge theory, respectively. 
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Figure 7: Average value of the real part of the plaquette, (Re Up), in the SU(144) TEK model 
versus the inverse bare 't Hooft coupling, h (from a hot start simulation). 
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Figure 8: Average value of the magnitude of traced link variables, (|;^Tr in the SU(144) 
TEK model versus the inverse bare 't Hooft couphng, b (from a hot start simulation) . 
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Figure 9: Typical profile of the eigenvalue densities of the link variables, U^, in the several 
Z^-breaking phases of the SU(144) TEK model; a labels the eigenvalues of SU(iV) matrices, 
which are pure phases of the form e*". 




Figure 10: Eigenvalue densities of a particular link variable at an inverse bare 't Hooft couphng 
oi b = 2.00, sitting at two different phases of the SU(144) TEK model, namely the twist-eater 
phase {Z'^^\ in the left) and the completely asymmetric phase {Z^, in the right). 
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Figure 11: Average value of the real part of the plaquette, (Re Up), in the original untwisted 
SU(81) Eguchi-Kawai model versus the inverse bare 't Hooft coupling, b. The nomenclature 
used for the different phases is the same as in [T3]. 
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